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P. ALBANO P. CANNARSA KHAI T. NGUYEN C. SINESTRARI 

Abstract. It is a generally shared opinion that significant information 
about the topology of a bounded domain Q, of a riemannian manifold M 
is encoded into the properties of the distance, dan, from the boundary 
of Q. To confirm such an idea we propose an approach based on the 
invariance of the singular set of the distance function with respect to 
the generalized gradient flow of of dan- As an application, we deduce 
that such a singular set has the same homotopy type as f2. 
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1. Introduction 

The distance function from the boundary of a bounded open set f2 C R n 
is a well-known object. Defined as 

ddn(x) = min \y — x\ Vx G R n , 

it is a nonsmooth function that occurs in different types of context. On 
account of its very definition, the distance retains some of the smoothness 
properties of the euclidean norm. For instance, dg^ is Lipschitz continuous 
in R™ with Lipschitz seminorm equal to one, and locally semiconcave in Q, 
that is, for every convex compact set G C !! there is a constant Kq £ R 
such that x h-> ddn(x) — Kg\x\ 2 /2 is concave on G. Moreover, the square 
of the distance is semiconcave with the uniform constant K = 2 in the 
whole space, that is, x i— > d 2 d ^{x) — \x\ 2 is concave in R n . Consequently, 
dgci is differentiable almost everywhere in R n . In particular, the distance is 
differentiable in f2 \ X, where X C SI is a set of Lebesgue measure zero called 
the singular set of don. In fact, general properties of semiconcave functions 
ensure that X is countably (n — l)-rectifiable, see [3]. 

The structure of the singular set has also been studied from the viewpoint 
of propagation of singularities. For general semiconcave functions, geometric 
conditions ensuring the propagation of singularities were obtained in [2]. 
Subsequently, in [3] and |BJ, singularity-propagation results were derived for 
viscosity solutions of the Hamilton-Jacobi equation 

F(x,u,Du)=0 in tt, (1.1) 

with F(x, u,p) convex in p. More precisely, given a singular point xo, which 
is not critical for u, one can show the existence of a nonconstant lipschitz 
arc 7 : [0, a[— > R n , starting at xq, which consists of points where u fails to 
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be differentiable. Moreover, 7 is a generalized characteristic of (jl.ip . that 
is, a solution of the differential inclusion 

Y(t) € co D p F( 7 (t), u( 7 (t)), £> + u( 7 («))) , (1.2) 

where D + -u denotes the superdifferential of u while 'co' stands for 'convex 
hull'. Observe that the above settings include the distance function, which 
solves the eikonal equation \Ddgn\ 2 = 1. In this case, the differential inclu- 
sion (|1.2p reduces to the gradient flow 7' € D + (1qq('j) — up to rescaling. 

The propagation of singularities along characteristics is a well-studied 
property of solutions to linear hyperbolic equations. In [10], for scalar hy- 
perbolic conservation laws in one space dimension, Dafermos observed that 
singular arcs could be regarded as generalized solutions of the same differ- 
ential equation governing the dynamic of classical characteristics. 

The above considerations can be naturally extended to an open subset f2 
of a riemannian manifold M. Here, in a local coordinate chart, the eikonal 
equation takes the form 

(A- 1 (x)Du(x),Du(x)) = 1 , (1.3) 

where A{x) is related to the riemannian scalar product g x on the tangent 
space T X M by the formula 

g x (Z,0 = (Mx)tO VC,CGT X M. 
In this case, the equation of generalized characteristics is 

7 , (t)G^- 1 ( 7 (t)) J D + n(7(0), (1-4) 

which will be referred to as the generalized gradient flow. 

The gradient flow of the distance function on a manifold has often been 
used in riemannian geometry as a tool for topological applications in con- 
nection with Toponogov's theorem, starting from the seminal paper [12] by 
Grove and Shiohama. A survey of the main results obtained by such tech- 
niques can be found in Chapter 11 of [IT]. However, these authors' approach 
differs from the one considered here, in that they used a regularization of 
the gradient flow of the distance which admits smooth solutions. 

In the aforementioned paper |10j . singularities are shown to propagate 
along a generalized characteristic, forward in time up to infinity. Intuitively 
speaking, such a behaviour is related to the well-known interpretation of the 
entropy condition for solutions of conservation laws (with a convex flux), 
which ensures that characteristics can only go inside a singularity: a char- 
acteristic which enters the singular set remains "trapped" there. 

Although the two-dimensional structure of the problem is essential for the 
proof of |1U| . one may wonder whether the same property holds for singular 
arcs of the distance function in arbitrary dimension. The results of this 
paper give a precise response to such a question showing that S is invariant 
for the generalized gradient flow, in the sense that any solution of (|1.4p . with 
7(0) G S, satisfies j(t) £ S for every t > 0. We note that such a property 
allows 7 to become constant at some time when a critical point of d has 
been reached. In order to better explain the main idea of our approach, we 
will treat the euclidean case first, in Theorem 13.71 below, and then extend 
the analysis to riemannian manifolds, in Theorem 14.51 
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The strategy of proof of our main results can be easily summarized, at an 
intuitive level, assuming that the arc 7 is smooth. For the euclidean distance 
in R™, basic results from [3] ensure that 7(5) G S if and only if |7'(s)| < 1. 
On the other hand, using properties of semiconcave functions we manage 
to show that the speed of a singular generalized characteristic satisfies the 
'logistic' differential inequality 



where K is a semiconcavity constant for Since the optimal constant 

for the distance function is K = 2, the above inequality yields |7'(s)| < 1 if 
the same is true for s = 0, forcing 7(5) G £. 

Surprisingly, the above method works in the riemannian case as well, even 
though dgQ fails to be semiconcave with K = 2, in general. One observation 
which is crucial for the extension is that, in fact, our idea is based on the 
existence of a nonlinear transform of the distance function which satisfies a 
suitable differential constraint. It turns out that cosh(a dg^), where a > 
is a suitable constant depending on the sectional curvatures of M, has the 
required properties, see Theorem 14.21 below. In addition, the proof requires 
the use of some standard tools in riemannian geometry, such as parallel 
transport and properties of geodesies. 

It is interesting to remark that semiconcave functions, and in particu- 
lar the distance function, have been studied in different domains, often 
independently. Applications of semiconcavity can be found in nonlinear 
partial differential equations ([12], [E]), riemannian manifolds and Alexan- 
drov spaces ([16]. [T8j). control theory ([7]), and optimal mass transporta- 
tion (|19|). Moreover, the singular set of the distance function is closely 
related to the cut-locus of the boundary of f2 — a widely studied object in 
riemannian geometry, mostly in the case of the distance from a single point 
(see e.g. [6]). 

More recently, interest in the properties of the singular set of dgn arose 
in applied domains such as computer science (see [15] and the references 
therein). Several authors have studied, at increasing levels of generality, the 
homotopy equivalence between S and f2. If the boundary of f2 is smooth (or 
piecewise smooth in dimension 2), then a homotopy can be constructed by 
moving every point outside S by the smooth gradient flow of the distance 
function until it reaches X. However, this method only works if the distance 
from a regular point p to £ along the gradient of the distance — the so called 
normal distance from £ — is a continuous function of p. If fl is a nonsmooth 
set of dimension at least 3, then this fails to be true and the construction of 
the homotopy becomes more involved [TS] • An essential step of the procedure 
is to extend the gradient flow of the distance function past a singular point, 
ensuring that the corresponding generalized characteristic stays singular for 
all time. This is why our result about the invariance of the singular set 
under such a flow can be used to provide an easy proof of the result of [15] 
and, more importantly, generalize the homotopy equivalence between £ and 
ft to the case of a complete riemannian manifold (see Theorem 15.31 below). 

As a possible application of the analysis developed in this paper, we would 
like to mention optimal exit time problems in R n that can be subsumed by 
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the riemannian setting. Specifically, we study a class of time optimal control 
problems for which the minimum time function T(x), which measures the 
minimum time needed to steer a point x E 0, to d£l, can be interpreted as 
a riemannian distance function (Example 15.41 below) . As a consequence, we 
deduce the homotopy equivalence between and the singular set of T(-), 
a result that would be hard to derive keeping the reasoning confined to 
euclidean space. 

This paper is organized as follows. First, in section [21 we introduce the 
essentials of our notation. Then, in section [3l we show the invariance of E 
in an open subset of M n , and generalize this result to riemannian settings in 
section 01 Finally, in section [5l we prove homotopy equivalence and discuss 
applications to optimal control problems. 



2. Notation 

Given x,y £ R n , we denote by (x,y) and by |ac| the euclidean scalar 
product and norm respectively. We set d(x,y) := \x — y\ for the distance 
between two points. 

If M is a riemannian manifold and x € M, we denote by T X M the tan- 
gent space to M at x and by T*M the cotangent space. There is a canonical 
isomorphism between the two spaces given by the scalar product. For sim- 
plicity, we use the same symbol (•, •) to denote the scalar product of two 
vectors of T X M, or of two elements of T*M, or the pairing of a form in T*M 
and a vector in T X M. The riemannian distance between two points x, y € M 
will be again denoted by d(x,y). 

For a function u defined in a subset of M n , we denote by Du and D 2 u 
the gradient and hessian of u where they exist. If u is a function on a 
riemannian manifold M, we denote by Du(x) £ T X M its gradient and by 
du(x) € T*M its differential; although the two notions are equivalent via the 
isomorphism recalled above, it will be convenient for our purposes to keep 
them distinct. In addition, we denote by D 2 u the hessian of u, interpreted 
as a linear operator from T X M to itself, as for example in [19} Ch. 14]. 

If u is a Lipschitz continuous function defined in a riemannian manifold, 
we denote by E(it) the set of points where u is not differentiable, and we 
call such a set the singular set of u. By Rademacher's theorem, E(it) has 
Lebesgue measure zero. 

If C is a nonempty closed subset of a riemannian manifold M (in par- 
ticular M could be W 1 ), we denote by dc(x) the distance function from C, 
defined as 

dc(x) = min d(y,x). (2-1) 

t/gC 

It is well known that dc(x) is Lipschitz continuous with constant 1 and that 
\Ddc{x)\ = 1 at every x where dc is differentiable. 

3. The euclidean case 

We first recall some properties of semiconcave functions. In this section, 
0, will be an open set in R n . 
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Definition 3.1. A function u : 0, — >■ R is called semiconcave if there exists 
K > such that 

I |2 

i-u(x) + (1 - i)u(y) - u(tx + (1 - i)y) < t(l - t)K ^ ~^ 

for any x, y G such that the segment from x to y is contained in O, and for 
any i G [0, 1]. We call K a semiconcavity constant for -u in O. We say that 
u is locally semiconcave in Q if it is semiconcave on any subset 4 CC 0. 

It is easy to see that u is semiconcave with constant K if and only if the 
function u(x) — -f-| x l 2 is concave or if D 2 u < K Id in the sense of distribu- 
tions, where Id denotes the identity matrix. 

The (Frechet) superdifferential of a function u : — > R at a point x G O 
is defined as the set 

^4- / \ f mn u(y) — u(x) — (p,y — x) 1 
J D+ U ( X ) = pel" : limsup H — <0 ■ 

In the case of a semiconcave function, the superdifferential enjoys the follow- 
ing properties; the proofs can be found in any textbook on convex analysis 
or in Chapter 3 of [7]. 

Proposition 3.2. Let u : f2 — > R be semiconcave. 

(i) The function u is locally Lipschitz continuous in £1 and differ 'entiable 
almost everywhere. 

(ii) The superdifferential D + u{x) is nonempty for all x G O. It is a 
singleton if and only if u is differ entiable at x, and in this case we 
have D + u(x) = {Du(x)}, where Du(x) is the standard gradient. 

(iii) For any x, y € such that the segment from x to y is contained in 
£1, for any p € D + u(x) and q G D + u(y), we have 

(q-p,y-x) <K\y-x\ 2 , (3.1) 

where K is a semiconcavity constant of u. 

(iv) Given {x n } C 0, such that x n — > x G Q and p n G D + u{x n ) such that 
Pn — > P, we have that p G D + u(x). 

The distance function dc provides an example of semiconcave function, 
as we recall here. 

Proposition 3.3. Given any nonempty closed set C C R n , the distance 
function dc is locally semiconcave on R n \ C . In addition, the squared 
distance function d c { ) is semiconcave on all R n with constant K = 2. 
Moreover, we have 

(d c {x)p-d c {y)q,x-y} <\x-y\ 2 (3.2) 
for all x,y G R n ,p G D + d(x) and q G D + d{y). 

Proof. The statements about the semiconcavity of dc and d c follow from an 
easy computation, see e.g. [7, Proposition 2.2.2]. Estimate ()3.2[) then follows 
from Proposition 13. 2f iii). observing that D + d c (x) = 2dc{x)D + dc{x). 

□ 
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In the following we will consider semiconcave functions u : Q — > R which 
solve an equation of the form 

(A- 1 (x)Du(x),Du(x)) = 1, ieH, (3.3) 

where A(x) a symmetric positive definite n x n matrix with C 1 dependence 
on x G S7. The formulation of the equation with A -1 , rather than A, is more 
convenient in view of the application to the distance function on riemannian 
manifolds in the next section, where A(-) will be the matrix associated with 
the metric on the tangent space. 

It is well known (see e.g. Prop. 5.3.1]) that, if u : — > M is semicon- 
cave, then the following properties are equivalent: 

• u satisfies (|3.3[) at every x £ Q where Du(x) exists; 

• for every x G f2 and p G D + u{x) we have (A~ 1 (x)p,p) < 1; 

• u is a viscosity solution of (|3.3[) (in the sense of [9]). 
Throughout the paper, we call a solution of (|3.3|) a locally semiconcave 
function satisfying the above properties. 

Given a solution of (|3.3p . we consider the differential inclusion 

7 '(t)G^- 1 ( 7 (t)) J D+ U ( 7 (t)). (3.4) 

A Lipschitz arc 7 : [0, to] — > is called a solution to the above problem 
if, for a.e. t G [0,ioL it satisfies 7'(i) = yl _1 (7(t))p(t) for some element 
p(t) G D + u(7(i)). Such an arc will also be called a generalized characteristic 
of equation (|3.3p associated with u. 

We now recall some properties of generalized characteristics. The main 
part of the statement (in particular claim (iv) about the propagation of 
singularities) follows from the results first proved in [3] and then obtained 
with a simpler approach in [8j 120] , For the convenience of the reader, we 
include the proof of some additional properties which were not explicitly 
observed in the above references. 

Theorem 3.4. Let u : f2 — > R be a solution of (|3.3p . Then, for every xq G f2 
there exists to > and a unique Lipschitz continuous arc 7 : [O^ot - ^ ^ 
which satisfies (|3.4p and the initial condition 7(0) = xq. In addition, the 
right derivative 7+(i) exists for every t G [0, to[, and p(t) := 7+(i) has the 
following properties: 

(i) p(t) G ^ 1 (7(t))L>+u(7(t)) /or every t G [0,t [ 
(p(*),A( 7 (t))p(t)) < (g,4(7(t))<7), Vg G yl- 1 ( 7 (t)) J D + n( 7 (t)). (3.5) 

(ii) p(i) is continuous from the right for every t G [0,io[ anc Z; f or oil 
points t* where it is discontinuous, we have 

Hmipf <p(t),A( 7 (t))p(t)) > (p(i*),A( 7 (t*))p(0). (3.6) 

(iii) For any t G [0,t [, 7 (i) G £(u) 2/ and on/y if (A(-f(t))p(t),p(t)) < 1. 

(iv) // xo G £(«), f/ien i/iere exists a G]0,to] smc/i t/iai 7 (t) G /or 
a// f G [0,a]. 

(v) For aZZ i G [0,to[ we have 

A.u( 7 (t)) = (A( 7 (t))p(t),p(t)), (3.7) 
where the symbol -4^ denotes the derivative from the right. 
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Before proving the theorem, we give an elementary continuous dependence 
result for generalized characteristics. 

Lemma 3.5. Under the above assumptions, for any U CC there exist 
C > and to > such that, if x,y G U and j x ,ly ar e solutions of t \3.4\ ) 
with initial conditions J x (0) = x and J y (0) = y respectively, then 

\lxit) - ly (t)\ < C\x - y\, te[0,t o }. (3.8) 

Proof. Let us take any ball B CC f2 and points x\,x 2 G B. Given any 
pi G D + u(x\) and p2 G D + u(x 2 ), we have 

(A~ 1 (xx)pi-A~ 1 (x2)p2,A(xi)(xi-X2)) = (pi-A(x 1 )A~ 1 (x 2 )p 2 ,x 1 -x 2 ) 

= (Pi -P2,a^i - £2} + {(I - A(xi)A(x 2 )~ 1 )p 2 , x\ - x 2 ) < Ci\x 1 - x 2 \ 2 , 

for some C\ > 0, where we have used property (hi) of Proposition 13.21 and 
the Lipschitz continuity of the map B 3 x i— > A~ 1 (x). 

Let now "f x ,ly be two solutions of (|3.4p contained in 5. Using the above 
estimate and the nondegeneracy of the matrix A(-) we hnd that 

^(A(lx(t))(lx(t) -7v(t)),7*(*) -7»(*)> 

^(7x(*))) (7*(t) -7v(<)),7*W -7y(*) 

+2( 7 ;(i) " 7^)), A(7*(t))(7» i (<) " 7,(0)) 

< C 2 \ lx (t)- ly (t)\ 2 

< C 3 (A( lx (t))(^ x (t) - ly (t)), lx (t) - 7y (t)). 
The Gronwall inequality yields 

(A( lx (t))( lx (t) - 7y(t)), lx (t) - 7y (t)) < C 4 \x - y\ 2 , 
which implies the conclusion. □ 

Proof, of Theorem \3.4\ The existence of an arc satisfying (|3.4p is proved in 
[HJ Theorem 3.2] (see also [3]), while uniqueness follows from Lemma 13.51 

The existence and continuity from the right of 7+(i) for every t, as well as 
properties (i) and (iv), follow from Corollaries 3.3 and 3.4 in [8]. Actually, 
those results require the additional assumption that ^ D + u{xq); however, 
if G D + u(xq) then the arc 7 is the constant one j(t) = xq, and all the 
properties of our statement are trivially satisfied. 

To prove (|3.6p . let us pick any sequence t n f t* such that lim n _ ) . 00 p(t n ) ex- 
ists. By (i), we have that A{^{t n ))p{t n ) G D + u( , y(t n )) for any n. Therefore, 
setting p_ = lim„_ s>00 p(i n ), we have by Proposition ^, 2f iv) that A(j(t*))p- G 
D + u{^f{t*)). But then we obtain from (|3,5p that 

(p_,A( 7 (t*))p_) > <p(t*),A( 7 (t*))p(t*)>, 
which implies (|3.6p and completes the proof of (ii). 

To prove (iii), suppose that u is differentiable at 7 (i). Then, by part (i) 
and Proposition I3.2f ii). we have that A{^{t))p{t) = Du{^{t)). Therefore, 
since u solves (|3.3p . (p{t), A(^(t))p{t)) = 1. Conversely, if u is not differen- 
tiable at 7(t), then D + u(j(t)) contains more than one point. As recalled 
previously, a solution u to equation (|3.3p satisfies {A{^{t))q^q) < 1 for all 
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q G D + u( , ~f(t)). Since D + u(7(t)) is a convex set containing more than one 
point and A is positive definite, we deduce that (A(^(t))q,q) < 1 for some 
q G D + u(j(t)). Then (ii) implies that (p(t), A(-y(t))p(t)) < 1. 

To prove property (v), we have to recall some details of the proof of the 
existence of the singular arc 7 given in [8j. The authors introduce there 
a family of smooth functions u k converging uniformly to u with bounded 
Lipschitz constant. The arc 7 is then obtained as the uniform limit of a 
sequence of smooth arcs "y k which solve the equation 

i k {t) = A-\ lk {t))Du k { lk {t)). 

The above properties easily imply that we also have weak convergence 7I — ! - 
7' in L 2 ([0, to], R n ). If we exploit the lower semicontinuity of convex func- 
tionals with respect to weak convergence, we obtain, for any < ti < t% < to, 

liminf f\A{ lk {t))i k {t),i k {t))dt 
k ^°° Jti 

= hminf /V(7(t))7^),7fc(*)Wi 
k ^°° Jti 

+ lim f\[A{ lk (t)) - A( 7 (t)M(t), 7 ' k (t))dt 

> / t2 (A(7(i))7 , (t),7 , (*))^ J 

where we have also used the uniform convergence of 7^ and the boundedness 
of 7I. It follows that 

u{^{t 2 )) - u{^{t{) = lim tifc(7fc(t 2 )) - Ukilkih)) 

= lim [ 2 (Du k ( 7k (t)), 7 ' k (t))dt= lim [ 2 (A( lk {t)) 1 ' k (t), 1 ' k (t))dt 

> r(A( 7 (t))7'(t),7'(t))dt = / <2 ^( 7 (t)Mt),p(t))dt, 
which implies 

um tt(7(t + ft) )-„ (7(t)) 

for every t € [0, to[, by the right continuity of p. On the other hand, since 
A(j(t))p(t) G D + u(j(t)), we also have 

uh(t + h)) - u(j(t)) , / . . . 7(t + h) - 7ft) \ 

limsup-^ -j KJ ±- LL < lim ( A h(t ))p(t , — r — ) 

h^o+ h h^o+ \ h I 

= (A(y(t))p(t),p(t)). 
This concludes the proof. □ 

Corollary 3.6. Let Q C 1" be an open set. Then, for every x G 0, there 
exists a unique Lipschitz continuous arc 7 : [0, oof— > £1 such that 

7 '(t) G D+d 9U (j(t)) t G [0, 00 [ a.e. 7 (0) = x. (3.9) 

Ln addition, for any to > such that 7(to) G S(dgn) t/iere exists a > suc/i 
t/iat 7(t) G S(cian) /or aZZ t G [to,io + Finally, the derivative from the 
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right 7+(i) exists for all t € [0, +00) and satisfies the properties described in 
Theorem\3.4[ with u(x) = dg^(x) and A(x) = Id. 



Proof. The statement follows directly from Theorem 13.41 provided we show 
that the maximal interval of existence of 7 is [0, +00 [. To see this we 
note that, if such an interval is [0, T[ with T 7^ +00, then necessarily j(t) 
approaches d£l as t — > T, that is, dgn(^(t)) — > as t — > T, in contrast with 
the property that docii^it)) is positive and nondecreasing in t by (|3.7p . □ 

We are now ready to prove the main result of this section. 

Theorem 3.7. Let C M n be an open set, let x € fl, and let 'j(-) be the 
solution of (|3.9p given by Corollary \3.b\ // 7(^0) £ £(<fefi) for some to > 0, 
then j(t) € T,(dgn) f or a ll t € [to, +00 [. 

Proof. For simplicity of notation, we suppose to = 0. Set 

p(t) := j' + (t), S(t) := 4fln(7(t)) t G [0,oo[. 

From Theorem 13.41 and Corollary 13.61 we know that p(t) € D + dgQ('y(t)) 
for all t; in addition j(t) € S(cign) an d b(*)l < 1 f° r all t > in a right 
neighbourhood of 0. Our aim is to show that |p(t)| < 1 holds for every t > 0; 
by part (iii) of Theorem 13.41 this will prove our assertion. 
Let < s < t. By Proposition 13.31 we have 

Ht)(p(t)-P(s),j(t)-J(s)) 

< | 7(t ) _ 7(s )|2 _ {6{t) _ «5( s ))<p( s ), 7 (t) - 7 ( s )). (3.10) 

Now, set t = s + h (h > 0), 

j(s + h)-j(s) 5(s + h)-5(s) 

j h (s) = and S h (s) = . 

n h 

We observe that I7/J < 1, |<5^| < 1, since 7 and 5 are both 1-Lipschitz 

function. From (|3.10p we obtain 

(Yftto, This)) < ^^[|T»i(*)| a -*hW<p(*),7hW>]. (3-11) 

Here and in the rest of the proof, we use for simplicity the notation of 
the ordinary derivative to mean the derivative from the right of expressions 
involving 7^. Observe that, by Theorem 13.41 we have for all s > 

lim n(s) =p(s), lim 5 h (s) = \p(s)\ 2 . (3.12) 

h— S-0+ ft— S-0+ 

From a heuristic point of view, it is useful to take the limit as h 4- in 
iKTTh . We obtain 

^|p( S )| 2 <^y|p( 5 )| 2 (l-b( S )| 2 ). (3.13) 

Such an inequality implies that, if |p(to)| < 1> then \p(s)\ < 1 for all s > to- 
However, such a reasoning is only formal, because we cannot say anything 
about the differentiability of p(-). It is interesting to observe that the crucial 
constant 1, in the expression (1 — |p(s)| 2 ) above, arises from the previous 
computations as K/2, where K = 2 is the semiconcavity constant of d^. 

Although the above argument is not rigorous, it suggests that 7h(-) can be 
estimated by a suitable adaptation of the separation of variables procedure 
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which could be used to integrate (|3. 13|) . To do this, let us first fix e > 
small. Since < |7/j(s)| < 1 for all s > 0, h > 0, we have that 

(| 7 h(s)| 2 + £)(! + £- hh(s)\ 2 ) >e(l + e) > 0, a > 0, h > 0. (3.14) 

Thus, we can divide both sides of (|3.1ip by the above expression to obtain 

(|7^)| 2 )' 2 \l h (s)\ 2 -6 h (s){p(s), lh (s)) 



(\ 7h (s)\* + e)(l + e- | 7/i ( S )| 2 ) " S(s + h) (\ lh (s)\* + e)(l + e - \ lh {s)\*) ' 
Integrating over [0,t], we find 



'" 3, .1+ E -|7,,(»)I 2 
where 



«(« + *) (h*MI 2 + e)(l+e-|7fcMI 2 ) 



Therefore, 



e^MW^ (i + e )_ e l + 2e 



where 

\lh(0)\ 2 + e 



-,h,e 



We now want to let first h | and then s I 0. By ()3.12p we find, 

So (fc5K ft,e ) ~ 40 1 + |p(0)| 2 ~ 1 - |p(0)| 2 ' 

On the other hand, since by ()3.14p the integrand in the definition of i*), e is 
uniformly bounded in h, we obtain, again by (|3.12p . 

kio v v *W (IpWI 2 + + e -bWI 2 ) 

< (l + 2e) / -^rds, 
Jo d(a) 

which implies 

limsup ( limFft £ (i) J = / — ^ (is =: a(t) . 

E+Q / JO <H S ) 

Thus, letting first /i J, and then e J, in (|3.15p . we conclude that 

m p < l^)! 2 ^ (3 16) 



Therefore, |p(i)| 2 < 1 for every i > 0, and the proof is complete. □ 
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4. The riemannian case 



Let us now consider a complete riemannian manifold M, possibly noncom- 
pact. To extend to this framework the techniques of the previous section, 
we first need to recall some basic properties of parallel transport and ge- 
odesic curves. On M, there is a canonical notion of derivative of a vector 
field, called covariant derivative. Using this definition, a vector field is called 
parallel along a curve 7 if its derivative in direction j'(t) is zero for all t. If 
we have a curve 7 : [a, b] — > M and a vector v G T 7 ( a )M, there is a unique 
vector field v(t) G XyU\M, with t G [a, b], which is parallel along 7; such 
a field v(t) is called the parallel transport of v along ~f(t). Parallel trans- 
port preserves the scalar product and therefore gives an isometry between 
the tangent spaces at different points. The geodesies on M can be defined 
equivalently as the curves 7 such that the speed j'(t) is parallel along the 
curve 7 itself or as the curves which are stationary for the energy functional. 
Geodesies have constant speed and are curves of minimal length between two 
endpoints if these points are close enough to each other. 

Given a point x G M, we denote by exp x (-) the exponential map at x. We 
recall that, given a tangent vector v G T X M, exp x (y) is the point reached at 
t = 1 by the geodesic 7(t) starting with 7(0) = x and 7'(0) = v. If / is a 
smooth function and df(x) G T*M is its differential at x, we have 

f(exp x (v))-f(x) = (df(x),v)+o(\v\), v£T x M, v -)• 0. 

Let us now consider a function u : M — > R not necessarily smooth. We 
say that p G T*M belongs to d + u(x), the superdifferential of u at x, if 

u(exp x (v)) — u(x) < (p, v) + o(|vj), v G T X M, v — > 0. 

This is equivalent to saying that there exists a smooth function / touching u 
from above at x such that df(x) = p. It is easy to see that, if p G d + u(x) and 
if 7 : [—a, a] — > M is any smooth curve such that 7(0) = x (not necessarily 
a geodesic), then 

limsup n(7W)-,(7(0)) ^ ^(0)). (4.1) 
/i->0 " 

We recall that a subset f7 C M is called convex if any distance mini- 
mizing geodesic between two points in U is contained in U. The notion of 
semiconcavity can be extended to riemannian manifolds as follows. 

Definition 4.1. A function u : U — > R, with U C M convex, is called 
semiconcave in J7 with constant if if, for every geodesic 7 : [0, 1] — > U and 
t G [0, 1], we have 

(1 - t)u( 7 (0)) + tu( 7 (l)) - u( 7 (t)) < *(1 - *) # d(7(0) 2 7(1))2 • (4.2) 

A detailed exposition of the basic properties of semiconcave functions 
on a manifold is given in [19]. Notice that, in such a reference, functions 
satisfying (|4.2p are called "semiconcave with modulus uj(t) = Kt 2 /2" . 

It can be checked (see Proposition 10.12 and inequality (10.14) in |19j) 
that, if u is semiconcave with constant K, then its superdifferential is 
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nonempty at each point. In addition, any p £ d + u(x) satisfies 




for all x £ U, v £ T X M such that exp x (v) £ U. 

Denote by 7(-) the geodesic j(t) = exp x (tv) starting at x with speed v. 
If we set y = 7(1) and w = 7'(1) £ T y M, we have x = exp y (—w). Thus, if 
q £ we find — < — (g, + if |u;| 2 /2. Since \w\ = d(x, y) = 

\v\, we can sum up with the previous inequality to obtain 



Denote by II : T X M — > T y M the parallel transport along the geodesic 7. 
Since parallel transport preserves the scalar product and 7 is a geodesic, 
we have w = Tlv and {q,w} = (II -1 (g), i;). We conclude that the above 
inequality can be rewritten as 



for all v £ T X M with ex.p x (v) £ U, and any p £ d + u(x) and q £ cT^i^exp^, (?;)). 

It is well known that the properties of the hessian of the distance func- 
tion in a riemannian manifold are closely related with the curvature of the 
manifold. Roughly speaking, positive curvature decreases the hessian of 
the distance function (i.e., gives a "stronger" semiconcavity) , while negative 
curvature increases it. In particular, it can be proved that the square of 
the distance function is semiconcave with constant 2 only if the manifold 
has nonnegative sectional curvature. Even in the case when the curvature 
has arbitrary sign, however, it turns out that we can replace the square by 
another function of the distance which enjoys the properties we need for our 
application. The crucial result for our purposes is the following. 

Theorem 4.2. Let M be a riemannian manifold and letVt C M be any open 
set (not necessarily smooth). Suppose that all sectional curvatures k at any 
point ofQ satisfy k > —a 2 for some a > and define v(x) = cosh.(adQn(x)) . 
Then, given any convex compact set C C £1, the function v is semiconcave 
on C with constant K = a 2 maxc v(x). 

Proof. In the case a = 1 the result is a direct consequence of Lemma 57 in 
Chapter 11 in [T7], where it is stated that, if the sectional curvatures of M 
are greater than —1 and if d p { ) is the distance function from a given point 
p £ M, then the function u(:r) = cosh(d p (x)) satisfies D 2 v(x) < v{x)Id 
everywhere. Since v(x) is not everywhere smooth, in general, the above 
bound on the hessian is understood in a suitable weak form [17} §9.3.1] 
which easily implies the stated semiconcavity estimate. Once the property 
is established for d p , it extends to the distance from an arbitrary set because 
the infimum of semiconcave functions with uniformly bounded semiconcavity 
constants is semiconcave with the same constant. The case of a general a is 
immediately reduced to this one by a homothety of the metric. Finally, we 
observe that the behaviour of the metric outside Q does not influence dgn(x) 
for x £ and therefore it suffices to assume the bound on the sectional 
curvature on Q. □ 




(4.3) 




(4.4) 
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Remark 4.3. We mention that, if the infimum of the sectional curvature 
is zero or positive, then it is possible to use functions different from the 
hyperbolic cosine which give a sharper semiconcavity estimate (e.g. the 
square of the distance in the euclidean case) ; the result of Theorem 14.21 
however, suffices for the purposes of this paper. This theorem also imples 
that the distance function itself is locally semiconcave in f2. However, the 
semiconcavity constant in general becomes unbounded as dO, is approached. 

We now proceed to show that Theorem 13.41 and Corollary 13.61 can be 
extended to manifolds. 

Theorem 4.4. For a given open bounded subset SI C M, let us set u(x) = 
dan{x) for x G 0. For every xo G O there exists a unique lipschitz continuous 
arc 7 : [0, +oo[— > Q, such that 

Y(t) G d + u{~f(t)) t G [0,+oo[ a.e. 7(0) = x . (4.5) 

The arc 7 satisfies properties analogous to the ones of Theorem \3.4\ and 
Corollary \ 3. 6\ in the euclidean case. In particular, the right derivative 7+(t) 
exists for every t > 0, is continuous from the right and satisfies (|4.5p every- 
where. The derivative of u along 7 satisfies 

^r«(7(*)) = |-/+(t)| a , *€[0,+oo[. (4.6) 

Moreover, for any to > such that 7 (to) S S(-u) there exists a = a (to) > 
such that j(t) £ ~S{u) for all t G [to, to + °~[ ■ 

Notice that, since 7'(t) E T 7(t) M and d + u(j(t)) C T* {t) M, in (gSJ the two 
spaces are identified via the canonical isomorphism. 

Proof. The result can be easily deduced from the euclidean case by using a 
local coordinate chart. In fact, if <f> : U — > M is a local chart around xq, 
where U C M n , and G(x) is the matrix associated to the scalar product on 
T X M in the chart cj>, then it is easy to see that the function u := uo(f> : U — > M 
satisfies 

1,3=1 J 

where g^(x) are the entries of the inverse matrix G _1 (x). Thus, the as- 
sertions of the theorem follow from the corresponding ones of Theorem 13.41 
We observe, in particular, that the equation satisfied by the generalized 
characteristics can be written in local coordinates as 

7 , (t)GG- 1 ( 7 (t)) J D + S(7(t)), (4.7) 

where D + u is the euclidean superdifferential of u. Finally, the property 
that 7(t) can be defined for t G [0, +00) is obtained by the same argument 
of Corollary □ 

We now show that Theorem 13.71 can be generalized to manifolds. 

Theorem 4.5. Let M be any smooth complete riemannian manifold, let 
f2 C M any bounded open set, let u(-) = dgn(-) and let j(-) be the arc of 
Theorem \4-4\ If j(po) £ ^( u ) f or some to > then j(t) G S(n) for all 
t G [to, +oo[. 
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Proof. Since Q is bounded, we can find a finite value a > such that the 
sectional curvature is everywhere greater than —a 2 on and Theorem 14.2 
can be applied. 

Let xo € ri and let 7 (i) be the solution of the differential inclusion (|4.5p . 
Let us set 

p(t)= 7 ' + (t), 6(t) = u( 7 (i)). 

We now suppose that 7(io) £ We then fix any s > io> and consider 

i > s close enough to s so that 7(s) and 7 (i) both belong to a neighborhood 
where any two points are connected by a unique minimal geodesic. 

Let us call v(s,t) the vector in T 7 ( S )M such that j(t) = exp 1 ^(v(s,t)). 
Also, we denote by Il Sit : T^r s \M — > T^mM the isometry induced by the 
parallel transport along the geodesic connecting 7(s) to j(t), and we set 
Tlt )S = H~l to denote the inverse map, which is associated to the same 
geodesic with the opposite direction. 

Using Gauss Lemma (see e.g. Lemma 3.3.5 in [H]), we obtain that 

d d( 7 ( S ),7(*)) 2 
dt+ 2 



(n S)t v( S ,i),7+(*)> = (v(s,t),U t>s p(t)), 



and, similarly, 



ds+ 2 
It follows, for h > small enough 

d d( 7 (s),7(» + '>)) 2 



8 <7W,7WJ! H _„ (M)i!)W) . 



ds+ 2 



n s+ftjS p(s + h) - p(s) , t>(s, s + h) ) . (4i 



Let us now set </>(r) = cosh(ar), for r£l. Then, Theorem 14.21 gives a 
semiconcavity estimate on the function eft o u. Let us also observe that, by 
the definition of superdifferential, 

q£d + u(x) <^=> <p'{u(x))q € d + {(p o u){x). (4.9) 

Let us denote by C C M the spherical neighbourhood of radius h centered 
at 7(5). If h is sufficiently small, C is a convex set which includes the point 
7(s + /i), and in addition sup c u(x) < u(7(s)) + /i. By Theorem l4.2l we have 
that (cj>ou)(x) is semiconcave on C with constant given by K = a 2 (j)(5(s)+h). 
Therefore, using (|4.4|) . (|4.5p . (|4.9p and the property that a 2 (j) = <fi", we find 

( <j>'(8(8 + h)) U s+htS p(s + h)- <f>'(6(s)) p(s) , v(s, s + h)} 
< ( j ) "(5(s) + h)d( 1 (s), 1 (s + h)) 2 . 

We rewrite the above inequality as 

4/(&{a + h)) (U s+hiS p{s + h)- p(s), v(s, s + h)) 
< 0"(5(s) + h)d(-f(s),~f(s + h)) 2 (4.10) 
-[4/(6(8 + h)) - 4/(6(8))] (p(s),v(s, s + h)). 

Let us set 

v(s,s + h) 6(s + h)-8(s) 
v h (s) = - h , 6 h (8) = . 
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Dividing inequality (|4,10p by h 2 , we obtain, using formula f|4. 8[) and the fact 
that d(7(s),7(s + h)) = \v(s, s + h)\ = h\v h (s)\, 

(4.11) 

<j/(6(s + h)) -<l/(6{ a )) 
h 

We rewrite for simplicity this inequality as 

d 



-(p(s),v h (s)). 



ds+ 



v h (s)\ 2 <Ms), (4.12) 



where we have set 



fi{6(s + h)) 

<j/(5(s + h))-<j/(6( a )) 
h 



J>"(6(s) + h)\v h (s)\ 2 

(p(s),v h (s)) 



(4.13) 



Now, we observe that 



hm^(s) = ~f' + (s) = p{s), (4.14) 



for a.e. s, which can be easily checked for instance by using local coordinates 
around 7(5). More formally, the above relation follows from the fact that 
linift^o Vh(s) = ^ \ h=0 v(s, s+h), that exp 7 ( s ) v(s, s + h) = ^(s + h) and that 
the differential of exp at zero is the identity (see the proof of Proposition 18 in 
Ch. 5 of [T7]); therefore, the derivatives J^| fe=0 v (s, s+h) and -^\ h=Q ^(s+h) 
coincide. 

In addition, using Theorem 14. 4( we find that for a.e. s 

lim <^ s + fr))-<//(^)) 
h^O h 

= ^( S ))lim / (S + ^- (5(s) = ^( S ))|p( g )| 2 . 

Therefore, we see that the function i^h{s) defined in (|4.13p is uniformly 
bounded for h > small and s varying in a bounded interval. Moreover, for 
all s > 0, 



(bWI 2 - W»)l 4 ) ■ (4-15) 



tanh(a 5(s)) 



From this point on, the proof proceeds as in the euclidean case. For e > 
sufficiently small, we obtain from (|4.12p 

(MS)|2) ' < , : Ma) 1 



(M»)l 2 + e)(l+e-K(»)l 2 ) " (h(.)P + £)(l+e-h( s )P) 
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for a.e. s > to- We assume for simplicity that to = 0. Integrating over [0,£], 
we have 

\v h (s)\ 2 +e n * 



log 



.1 + e - \v h (s)\ 2 
where 

iph(s 



q < */,, e (t) 



* fc , E (t) = (l + 2e) / 
•/ o 

Therefore, 



(K( S )|2 + £ )(l + e -|^( s )| 



2^ 



ds . 



< : - T h a , 7 ■ (4-17) 

1 + e ^h,sW^ h 



.,£ 



where 



c . k(0)| 2 + e 

C,h,e " 



l + e-M0)| 2 ' 
Using (|4.15p we obtain 

lim ( lim $/j ^] = f 2 " .. ds =: (3(t). 

eio \hio ' ) J tanh(ao(s)) 

Letting first h I and then e 1 in (|4.17p we obtain, thanks to (|4.14p , 

" |p(0)| 2 e^W + l-|p(0)| a ' 
which implies that |/>(t)| 2 < 1 for every t > 0. □ 

5. HOMOTOPY EQUIVALENCE 

Let us begin by recalling the well-known notion of homotopy equivalence. 

Definition 5.1. Let X and Y be two topological spaces and let 

f : X ^Y and g : X ->• Y 

be two continuous maps. We say that / and <? are homotopic if there exists 
a continuous map iJ : X X [0, 1] — > Y, called homotopy, such that 

H(0, •) = /(■) and £T(1, •) = <?(•). 

Furthermore, we say that X and Y have the same homotopy type if there 
exist continuous maps / : X — > Y and g : Y — > X such that jo / and / oj 
are homotopic to the identity on X and Y, respectively. 

The following result is a direct consequence of Definition 15.11 

Lemma 5.2. Zet Y (Z X. If there exists a continuous map 

H : X x [0, 1] -»■ X 

suc/i t/iat 

(a) -ff (x, 0) = x, /or every x € X, 

(b) i?(x, 1) € Y , for every x G X, and 

(c) ff(x, t) G Y ; /or every (x, t) G Y x [0, 1], 
t/ien X and Y have the same homotopy type. 

We are now in a position to apply Theorem 14.51 in order to obtain the 
following homotopy equivalence result. 
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Theorem 5.3. Let f2 be a bounded open subset of a smooth riemannian 
manifold M. Then Vt has the same homotopy type as E(dgn). 

Notice that the above theorem requires no regularity assumption on <9$7. 

Proof. In view of Lemma l5?2l it suffices to construct a continuous map 

H : Q x [0, 1] -> O 

satisfying conditions (a), (6), and (c) above with Y = E(dan). 

For any x G f2, let 7a;(-) be the generalized characteristic starting at x. 
We claim that 

3T>0 : Vxefi 7ie (T) G E(tfen)- (5-1) 
Indeed, set T = 2diam(f2) and let x G Q. Arguing by contradiction, suppose 
j x (T) i S(d an ). Then, in light of Theorem WM lS) £ Z(d dn ) for all 
t € [0,T]. So, |7^.(i)| = 1 for every t G [0, T\. Hence, owing to (|4.6p . we have 

dan(lx(T)) = d dn (x) + / \^' x (t)\ 2 dt = d 9n (x) + T. 

Jo 

Then, 

2diam(0) = T = ddn(j x (T)) - don(x) < diam(fi). 
The above contradiction shows that (|5.ip holds true. Next, define 
H{x,t) = lx {tT) (i,t)efix[0,l], 

We point out that H is a locally Lipschitz continuous map in view of 
Lemma 13.51 Moreover, on account of (|5.1|) and Theorem 14.5} H satisfies 
conditions (6) and (c) of Lemma 15.21 This completes the proof. □ 

The analysis of this paper on riemannian manifolds applies, in particular, 
to optimal exit time problems in R n . In the example below, we study a 
time optimal control problem, deducing a homotopy equivalence result that 
would be hard to obtain arguing just inside the euclidean framework. 

Example 5.4. Let F : W 1 — > M nxn be a smooth function such that det F(x) ^ 
for all x and let £1 C M. n be a bounded open set. For any given x £ £1 we 
consider the control system 

( y'(t) = F(y(t))a(t), 

1 y(o) = x, 

where a : [0, +oo) — > -Bi(O) is a measurable function called the control. We 
denote by y(-;x, a) the trajectory of (j5.2|) . and we define the exit time from 
of the trajectory as 

r(x,a) = inf{i > : y(t;x, a) G dO,} G (0, +oo]. 

The minimum time function is defined as 

T(x) = inf t(x, a), i£fi. 

a 

Under our hypotheses, it is well known that the infimum is attained and that 
T(-) is a semiconcave solution of the Hamilton-Jacobi-Bellman equation 

H{x,DT(x)) = 1, i£fi, (5.3) 

where H is defined as 

H(x,p) = (F(x)F*(x)p,p) 
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with F* the transpose matrix. 

Let us consider the riemannian metric g on W 1 induced by the scalar 
product with matrix 

G(x) := (F*)~ 1 (x)F~ 1 (x). 

Then, using the subscripts e and g to distinguish between the euclidean and 
riemannian metrics, we have 

\v\ g <l (G{x)v,v) e < 1 ^(x^U < 1, 

which shows that an arc y(-) is an admissible trajectory for the control 
system (|5.2p if and only \y'(t)\ g < 1. It follows that T(x) = dgn(x), where 
the distance function dg^i is taken with respect to the riemannian metric g. 

Thus, the previous analysis can be applied to the singular set £(T) of the 
minimum time function. In particular, recalling also (|4.7p . we obtain that 
£(T) is invariant under the flow induced by the differential inclusion 

7 '(i) € G- 1 ( 7 (t))D+T( 7 (t)). 

The above inclusion, up to a factor 2, can be written equivalently as 

7 / (t)€Z> p ff( 7 (t),D+r(7(i))), 

which is the equation of the characteristics associated with (|5.3|) . Moreover, 
Theorem 15.31 ensures that and S(T) have the same homotopy type. 
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